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Abstract 

In this paper, we consider elliptic differential operators on compact 
manifolds with a random perturbation in the 0th order term and show 
under fairly weak additional assumptions that the large eigenvalues 
almost surely distribute according to the Weyl law, well-known in the 
self- adjoint case. 

Resume 

Dans ce travail, nos considerons des operateurs differentiels ellip- 
tiques sur des varietes compactes avec une perturbation aleatoire dans 
le terme d'orde 0. Sous des hypotheses supplementaires assez faibles, 
nous montrons que les grandes valeurs propres se distribuent selon la 
loi de Weyl, bien connue dans le cas auto-adjoint. 
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1 Introduction 



This work is a continuation of a series of works concerning the asymptotic 
distribution of eigenvalues for non-self-adjoint (pseudo-) differential opera- 
tors with random perturbations. Since the works of L.N. Trefethen [8], 
E.B. Davies [2], M. Zworski [9J and many others (see for instance jl] for 
further references) we know that the resolvents of such operators tend to 
have very large norms when the spectral parameter is in the range of the 
symbol, and consequently, the eigenvalues are unstable under small pertur- 
bations of the operator. It is therefore quite natural to study the effect of 
random perturbations. Mildred Hager [1] studied quite general classes of 
non-self-adjoint /i-pseudodifferential operators on the real line with a suit- 
able random potential added, and she showed that the eigenvalues distribute 
according to the natural Weyl law with a probability very close to 1 in the 
semi-classical limit (/i — > 0). Due to the method, this result was restricted 
to the interior of the range of the leading symbol p of the operator and with 
a no n- vanishing assumption on the Poisson bracket {p, p} . 

In [5] the results were generalized to higher dimension and the boundary 
of the range of p could be included, but the perturbations where no more 
multiplicative. In [6l [7j further improvements of the method were introduced 
and the case of multiplicative perturbations was handled in all dimensions. 

W. Bordeaux Montrieux [T] studied elliptic systems of differential oper- 
ators on with random perturbations of the coefficients, and under some 
additional assumptions, he showed that the large eigenvalues obey the Weyl 
law almost surely. His analysis was based on a reduction to the semi-classical 
case (using essentially the Borel-Cantelli lemma), where he could use and ex- 
tend the methods of Hager [3j. 

The purpose of the present work is to extend the results of [T] to the case 
of elliptic operators on compact manifolds by replacing the one dimensional 
semi-classical techniques by the more recent result of [7]. For simplicity, we 
treat only the scalar case and the random perturbation is a potential. 

Let X be a smooth compact manifold of dimension n. Let P° be an 
elliptic differential operator on X of order m > 2 with smooth coefficients 
and with principal symbol p(x, ^). In local coordinates we get, using standard 
multi-index notation, 

p'=J2 «°(^)^"' p(^'^)= E «°(^)^" (1-1) 

|Q!|<m |a|=m 

Recall that the ellipticity of P*^ means that p{x,C,) 7^ for ^ 7^ 0. We assume 
that 

p{T*X) ^ C. (1.2) 
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Fix a strictly positive smooth density of integration dx on X, so that the 
norm || ■ || and inner product (-j ■ ■) are unambiguously defined. Let 
r : L?'{X) L'^{X) be the antilinear operator of complex conjugation, given 
by Tu = u. We need the symmetry assumption 

P* = rPr, (1.3) 

where P* is the formal complex adjoint of P. As in [7J we observe that the 
property (11.31) implies that 

p(x,-0 =p(a;,0, (1-4) 

and conversely, if (11.41) holds, then the operator |(P + FPF) has the same 
principal symbol p and satisfies (II. 3p . 

Let R be an elliptic differential operator on X with smooth coefficients, 
which is self-adjoint and strictly positive. Let eo,ei,... be an orthonormal 
basis of eigenfunctions of R so that 

Re, = {^i]ft,, < /i° < /i? < < ... (1.5) 

Our randomly perturbed operator is 

P° = P + g°(a:), (1.6) 

where uj is the random parameter and 

oo 


Here we assume that Q;j(u;) are independent complex Gaussian random vari- 
ables of variance and mean value 0: 

a°~Ar(0,cr|), (1.8) 

where 

{^i])-'e-^^^^^' < a, < (;.;)-^ (1.9) 

3n — - 1 
M = , < /3 < -, p > n, (1.10) 

s ^ t / 



where s, p, e are fixed constants such that 



n n n 

- < s < p , 0<e<s . 
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Let H^{X) be the standard Sobolev space of order s. As will follow from 
considerations below, we have g° G H'^^X) almost surely since s < p — f • 
Hence g° G L°° almost surely, implying that P° has purely discrete spectrum. 

Consider the function F{uj) = a.Tgp{uj) on S*X. For given 6'o G ~ 
R/(27rZ), A^o e N := N \ {0}, we introduce the property P(^o, Nq): 

No 

J2 IV'FHI ^ on {cj G S*X- F{uj) = 9o}. (1.11) 

1 

Notice that if P{6q,Nq) holds, then P{6,Nq) holds for all 6 in some neigh- 
borhood of 9q. 

We can now state our main result. 

Theorem 1.1 Assume that m > 2. Let Q < 9i < 62 < 2tt and assume that 
P^Oi, No) and P(^2, A'o) hold for some Nq G N. Let g e C°°{[ei, 62]; ]0, oo[) 
and put 

n,,er,o,x = {^e^'; 0^<9<92, 0<r< XgiO)}. 

Then for every 6 g]0, | — /3[ there exists C > such that almost surely: 
3C{uj) < 00 such that for all X G [1, oo[.- 

nr^,,,o,A) - (^^oip"'n,,,o,A)l (1-12) 

Here clP^) denotes the spectrum and #(v4) denotes the number of elements 
in the set A. In 1^1.12) the eigenvalues are counted with their algebraic mul- 
tiplicity. 

The proof actually allows to have almost surely a simultaneous conclusion 
for a whole family of 6i,62,g'- 

Theorem 1.2 Assume that m > 2. Let Q be a compact subset of [0,27r]. 
Let No ^ N and assume that P{6,No) holds uniformly for 9 E Q. Let Q be 
a subset of {{g, 61,62); 6j G 0,0i < 62, g G ^2]; ]0, cx)[)} with the 

property that g and 1/g are uniformly bounded in C°°([^i, ^2]; ]0, oo[) when 
(g, ^1; ^2) varies in Q . Then for every S g]0, |— /?[ there exists C > such that 
almost surely: 3C{uj) < 00 such that for all AG [1, C)o[ and all {g, 61, 62) G Q, 
we have the estimate lil.lB) . 

The condition (11. 9p allows us to choose aj decaying faster than any neg- 
ative power of fi^j. Then from the discussion below, it will follow that quj{x) 
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is almost surely a smooth function. A rough and somewhat intuitive inter- 
pretation of Theorem 11.21 is then that for almost every elliptic operator of 
order > 2 with smooth coefficients on a compact manifold which satisfies the 
conditions (11. 2p . (11.31) . the large eigenvalues distribute according to Weyl's 
law in sectors with limiting directions that satisfy a weak non-degeneracy 
condition. 

2 Volume considerations 

In the next section we shall perform a reduction to a semi-classical situation 
and work with h^^Po which has the semi-classical principal symbol p in (11.11) . 
As in |5l El [7j , we introduce 

V,{t) = vol {p e T*X; \p{p) - z\'^ <t}, t> 0. (2.1) 

Proposition 2.1 For any compact set K G C = C \ {0}, we have 

V,{t) = 0(r), uniformly for z e K, < t < 1, (2.2) 

with K = 1/2. 

The property (12. 2p for some k g]0, 1[ is required in [5l El [7j near the 
boundary of the set F, where we count the eigenvalues. Another important 
quantity appearing there was 

vol(7 + /^(0,t)), (2.3) 

where ^ = dT and F d C is assumed to have piecewise smooth boundary. 
From (12. 2p with general k, it follows that the volume (12. 3p is 0{t'^'^~^), which 
is of interest when k > 1/2. In our case, we shall therefore investigate 
vol (7 + -B(0,t)) more directly, when 7 is (the image of) a smooth curve. The 
following result implies Proposition 12.11 

Proposition 2.2 Let 7 be the curve {re*^ G C; r = g{0), 9 G S^}, where 
Q<ge C\S^). Then 

vol(p-i(7 + /^(0,t))) = O(t), t^O. 

Proof. This follows from the fact that the radial derivative of p is 7^ 0. More 
precisely, write T*X \ 3 p = ru, u G S*X, r > 0, so that p{p) = r"^p{u), 
p{uj) 7^ 0. If p G p^^{j + D{0, t)), we have for some C > 1, independent of t, 

g{a.Tgp{uj)) -Ct< r'"|p(^)| < ^(argp(^)) + Ct, 
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/ g{aigp{uj)) -Ct \ - / g{aigp{uj)) + Ct \ 

so for every u G S*X, r has to belong to an interval of length 0{t). □ 

We next study the volume in fl2.3p when 7 is a radial segment of the form 
[ri, r2]e*^°, where < ri < r2 and Oq E S^. 

Proposition 2.3 Let 9q E , Nq E ISi and assume that P{9q,Nq) holds. 
Then if < ri < r2 and 7 is the radial segment [ri,r2]e'^^° , we have 

vol(p-i(7 + /^(0,t))) = O(ti/^«), t^O. 

Proof. We first observe that it suffices to show that 

yo\s*xF-\[9o -t,9o + t]) = 

This in turn follows for instance from the Malgrange preparation theorem: 
At every point uq E F~^{9o) we can choose coordinates ui, co'2n-i, centered 
at uo, such that for some k E {1, A'o}, we have that dl_^{F — 9o){uo) is = 
when < j < A; — 1 and 7^ when j = k. Then by Malgrange 's preparation 
theorem, we have 

F{uj) -9o = G{UJ){UJ^ + ai{uJ2, ...,UJ2n^l)(^l~^ + ••• + ak{uJ2, ■■■,UJ2n-l)), 

where G, aj are real and smooth, G{ujq) 7^ 0, and it follows that 

vol {F~\[9o -t,9o + t])n neigh M) = 0{t'/'). 

It then suffices to use a simple compactness argument. □ 
Now, let <9i < 92 <27c, g E C°°{[9i, ^2]; ]0, oo[) and put 

^Ler,rur2 = {^e^'' Oi<0< O2, ng{9) < r < r2^?(^)}, (2.4) 

for < ri < r2 < 00. If < ri < r2 < +cx) and P{9j, Nq) hold for j = 1, 2, 
then the last two propositions imply that 

yo\p-\dTl,^.^,^^^^ + D(0,t)) = O(tV^o)^ t ^ 0. (2.5) 

3 Semiclassical reduction 

We are interested in the distribution of large eigenvalues C of P°, so we make 
a standard reduction to a semi-classical problem by letting < h '^1 satisfy 

C=T^, /.xlCr^/-, (3.1) 
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and write 



h'^iP!^ - C) = h'^P!^ -z=:P + K^ql - z, (3.2) 

where 

p = hJ^p^ = J2 aaix; h){hDY. (3.3) 

\a\<ra 

Here 

a„(x;/i) = C(/i™-l"l) in (3.4) 
aa{x;h) = a°(x) when |a| = m. 

So P is a standard semi-classical differential operator with semi-classical 
principal symbol p{x,^). 

Our strategy will be to decompose the random perturbation 

where the two terms are independent, and with probability very close to 1, 
SQuj will be a semi-classical random perturbation as in [7] while 

\\ku;\\H'' < h, (3.5) 

and 

n, , . 

s^]^,P-^[ (3.6) 

is fixed. Then /i™P° will be viewed as a random perturbation of /i™P°+/C(^. In 
order to achieve this without extra assumptions on the order m, we will also 
have to represent some of our eigenvalues a°(co') as sums of two independent 
Gaussian random variables. 
We start by examining when 

Wh-'q'jH^ < h. (3.7) 



Proposition 3.1 There is a constant C > such that ( fj*. ly holds with prob- 
ability 

Proof. We have 

oo 
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3J 1^ 



(3.8) 



and the aj are independent. Now, using standard functional calculus for R 
as in [HI E] , we see that 

oo 

ll^'"?°lll«^El(/^?)^«^MI'' (3.9) 



where {n^Yaj ~ A/'(0, are independent random variables and = 

Now recall the following fact, established by Bordeaux Montrieux [T], 
improving and simplifying a similar result in [5]: Let c?o,(ii,... be a finite 
or infinite family of independent complex Gaussian random variables, dj ~ 
A/'(0, < 5j < CX3, and assume that < oo. Then for every t > 0, 

nY. > ^) < -p(^^^(^ - ^0 E^J))- (3-10) 

Here P(v4) denotes the probability of the event A and Co > is a universal 
constant. The estimate is interesting only when t > CoY^a^ and for such 
values of t it improves if we replace {do, di, ...} by a subfamily. Indeed, J2 ^"j 
will then decrease and so will maxa|. 

Apply this to (13. 9p with dj = {jJ^Yaj, t = h? . Here, we recall that 
aj = {^yK^(Tj, and get from ([II]), that 



max?| X /i^"^, (3.11) 



while 



oo 







Let iV(/i) = #(cr(v i?)n]0,/x]) be the number of eigenvalues of V in ]0,yu], 
so that iV(/i) X /i"' by the standard Weyl asymptotics for positive elliptic 
operators on compact manifolds. The last sum in (13.121) is equal to 



fi^^'-p'^dN{fi) = 2(p - s)fi^^'-^^-^N{fi)dfi, 
Jo 

which is finite since 2{s — p) + n < hj (13.61) . Thus 

oo 


and the proposition follows from applying (13.91) . (13. lip . (I3.12p to (I3.10p with 
t = h'^. □ 
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We next review the choice of parameters for the random perturbation in 
[7] (and [U]). This perturbation is of the form 5Q^^, 

= h^'^q^, 6 = ro/i^^+", < ro < v^, (3.14) 

where 

Qu^ix) = ^ ak{uj)ek{x), \a\cD < R, (3.15) 

0<hfj.l<L 

and a possible choice of L, i? is 

L = Ch-^, R = Ch-^, (3.16) 

with 

3n — K ~ 3n ,n , „x 

M = — M = — -K + {- + e)M. (3.17) 

Here e > is any fixed parameter in ]0, s — ^[ and k, g]0, 1] is the geometric 
exponent appearing in (12.21) . in our case equal to 1/2. 
The exponent A^i is given by 

' — ' 71 

Ni = M + sM + -, (3.18) 

and should be subject to a probability density on Bqd{0, R) of the form 

C(/i)e*(";'^)L(da), where 



Va$| = 0{h-^'), (3.19) 



for some constant A^4 > 0. 
Write 



q^ = ql + ql (3.20) 
= E «?(^)^^' = E ^»'r (3.21) 

From Proposition 13.11 and its proof, especially the observation after (13.101) . 
we know that 

Wh'^q'jHs < h with probability > 1 - exp(Co - ^^J^(^)- (3-22) 
We write 

P + h"^q'^ = iP + h'^ql) + h"^ql 

and recall that the main result in jTj is valid also when P is replaced by the 
perturbation P + h'^q'^, provided that < h. 
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The next question is then wether /i^g^ can be written as r^h'^^^^^q^ where 
q^^ = ^o<hAtO<L "^i^i lo^lc^ ^ with probabihty close to 1. We get 



1 

^0 



^0 To 

Applying (13.101) . we get 

P(l«lc. > i?^) < exp(C - ^^,(f_;V„) ), (3.23) 
which is 0(1) exp(— /i^'^) provided that 

-2M + 2 ^'^^}Jj'^} + 2{2N, +n-m)<-6. (3.24) 
m[l/nj 

Here tq < Vh and if we choose tq = or more generally bounded from 
below by some power of h, we see that (I3.24p holds for any fixed 6, provided 
that m is sufficiently large. 

In order to avoid such an extra assumption, we shall now represent a° 
for hfi^ < L as the sum of two independent Gaussian random variables. Let 
Jo = jo(^) be the largest j for which hjjP- < L. Put 

a' = ^h^e-^^~\ where K > p{M +1), C> 1 (3.25) 

so that a' < \aj for 1 < j < jo{h). The factor is needed only when 
(3 = 0. 

For j < jo, we may assume that a°(co') = cn'j^uj) + a'jiuj), where a'j ~ 
A^(0, (o"')^), a'j ~ A/'(0, {(Tj)^) are independent random variables and 

(^j = ) + {(^j > , 

so that 

< = - (^')^ X a,. 
Put ql = qi + q'^, where 
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Now (cf we write 

p + h^ql = {P + h^{q': + ql)) + h^q'^. 

The main result of |7| is valid for random perturbations of 

P,:=P + h"^{q: + ql), 

provided that \\h"^{q" + q'^)\\H^ < h, which again holds with a probability 
as in (13.221) . The new random perturbation is now h'^q'^ which we write as 

Toh'^'^'^~^"-q^, where g^j takes the form 

U^)= Yl /^^H^^' (3-26) 
with new independent random variables 

To ■' To 

Now, by fl3J0D . 



Vmlo > R") < exp(0(l)D ^'""o 



Here by Weyl's law for the distribution of eigenvalues of elliptic self-adjoint 
differential operators, we have D x [L/h)^. Moreover, L,R behave like 
certain powers of h. 

• In the case when /3 = 0, we choose tq = h^/'^. Then for any a > we 
get 

CD > -R) < Cexp( 



for any given fixed a, provided we choose K large enough in (13.251) . 

• In the case /3 > we get the same conclusion with tq = h~^a' if K is 
large enough. 

In both cases, we see that the independent random variables [3j in (13.261) . 
(I3.27P have a joint probability density C{h)e'^^"''^^ L{da), satisfying (13.191) for 
some depending on K. 

With K = 1/2, we put 

e,{h) = h^i{\n]-y + \n-), 

h To 

where tq is chosen as above. Notice that eo{h) is of the order of magnitude 
^«;-/3 ^ power of In ^. Then Theorem 1.1 in [7J gives: 
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Proposition 3.2 There exists a constant A^4 > depending on p,n,m such 
that the following holds: Let F d C have piecewise smooth boundary. Then 
3C > such that for 0<r<l/C, ?> Ceo{h), we have with probability 

^ ^ CeoW _g-c^ _ ce-d,, (3.28) 

^^n+max(n(M+l),iV4+M) ^ ' 

that 

|#(/.™p°)nr)-^^voi(p-i(r))| < (3.29) 

^(1 + C{r + ln(-)vol {p'\dr + D{0, r))))). 
n"- r r 

As noted in [Hj this gives Weyl asymptotics provided tliat 

ln(-)volp~^(9r + D(0, r)) = C(r"), (3.30) 
r 

for some a €]0, 1] (wliicli would automatically be the case if k had been larger 
than 1/2 instead of being equal to 1/2), and we can then choose r = e^/(i+"), 
so that the right hand side of (13.291) becomes < Cei+^/i"". 

As in [HI [7j we also observe that if F belongs to a family Q of domains 
satisfying the assumptions of the Proposition uniformly, then with probabil- 
ity 

>1 Ceo{h) 

^2^n+max(n(Af+l),Af4+M) 

the estimate (I3.29P holds uniformly and simultaneously for all F G ^. 



4 End of the proof 

Let 61,62, No be as in Theorem [HH so that P{6i,No) and P{62,Nq) hold. 
Combining the propositions 12.11 12.21 12.31 we see that (I3.30p holds for every 
a < I/Nq when F = F^^ q^.^ ^, A > fixed, and Proposition 13.21 gives: 

Proposition 4.1 With the parameters as in Proposition \3.2\ and for every 
a e]0, -^l, we have with probability 

> 1 _ _ ^M^) ^e-^ - Ce-^ (4.1) 

^1+^ /^ra+max(n(M+l),Af4+Af) 

that 

\#Hh-P^)r^n.,er,i,x) - (^^ol(p-^n,..;i,A))l < C'-^- (4.2) 
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Moreover, the conclusion ^4-2^ is valid simultaneously for all A G [1, 2] and 
all {9i, 62) in a set where P{9i, Nq), P{92, Nq) hold uniformly, with probability 
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>1 ^^^^ ^e"^-Ce-^. (4.3) 

For < 5 < 1, choose e = h-^eo < Ch^-^~\\nlf, so that e/eo = h'^. 
Then for some N^, we have for every a &]0,1/Nq[ that 

(4.4) 

simuhaneously for 1 < A < 2 and all {01,62) in a set where P{6i,Nq), 
P{62,No) hold uniformly, with probability 

>l-^e-^. (4.5) 

Here a/{l + a) / l/(A^o + 1) when a / 1/Nq, so the upper bound in fjO|) 
can be replaced by 

^^(i-/3-25)/(JVo+l) 

Assuming P[9i, Nq), P{62, Nq), we want to count the number of eigenval- 
ues of P^ in 

r — v3 

1,A — -L 6li,e2;l,A 

when A 00. Let k{X) be he largest integer k for which 2^ < A and 
decompose 

A:(A)-1 

Tl.A = ( [J r2fc,2'= + U r2fc(A) 


In order to count the eigenvalues of P° in r2fc_2'=+i we define h by = 1, 
h = 2"'=/™, so that 

#(a(p°) n r2.,2^+0 = #(a(/i'"P°)nri,2), 

-VOl(p"^(r2fc,2*+i)) = — — — -voi(p"^(ri,2)). 



(27r)" ' ^ (27r/i)" 

Thus, with probability > 1 — C2' m e"^™/^ we have 

|#(a(p°)nr2.,2^+0-7TA^voij,-i(r2. ,2^+01 < c,2^2-^(^-^-^^)^. (4.6) 
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Similarly, with probability > 1 - C2^skix)/m^-2^'^w/^"/c ^ ^lave 

mcr{P'Jnr,,^.^~,)-j^voip-\r,,^^^^^^^^^ (4.7) 

simultaneously for all A G [A,2A[. 

Now, we proceed as in [1], using essentially the Borel-Cantelli lemma. 
Use that 

oo 

e 

to conclude that with probability > 1 — C2'^5^e"^*'"/'^, we have 

|#(a(p°) n r2.,A)l < CsX"^'"^''^^^ + c{u) 

for all A > 2^. This statement implies Theorem ll.il □ 
Proof of Theorem 11.21 This is just a minor modification of the proof of 
Theorem 11.11 Indeed, we already used the second part of Proposition 13.21 to 
get fl4.7p with the probability indicated there. In that estimate we are free 
to vary {g, 9 1,62) in Q and the same holds for the estimate (14. 6p . With these 
modifications, the same proof gives Theorem 11.21 □ 
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